The propagation of Elastic Waves in Granular Solid Hydrodynamics 
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The anisotropic, stress-dependent velocity of elastic waves in glass beads - as observed by Y. Khi- 
das and X. Jia, see [Phys. Rev. E, 81:021303, Feb. 2010] - is shown to be well accounted for by 
"granular solid hydrodynamics," a broad-range macroscopic theory of granular behavior. As the 
theory makes no reference to fabric anisotropy, the influence of which on sound is in doubt. 

PACS numbers: 91.60.Lj, 81.40.Jj 



I. INTRODUCTION 

The anisotropic spectrum of elastic waves in granu- 
lar media is studied in [1]. Wave velocity was measured 
under varying stress, in glass bead samples that have un- 
dergone two different preparations: Rain deposition that 
consists of pouring the glass beads into the acoustic cell 
through a mesh; and de- compaction, implying the gentle 
motion of a horizontal grid through the bead pack, from 
the bottom to the top, after the deposition. As the veloc- 
ity was seen to depend on the preparation, and because 
it is known from numerical simulations that the rain de- 
position creates an anisotropic distribution of the contact 
angles, with two preferred directions oriented roughly at 
30° from the vertical, this was taken as evidence of the 
influence of fabric anisotropy. 

However, also reported was the fact that the prepa- 
rations lead to different densities. So a crucial question 
is whether the density alone suffices to explain the ob- 
served difference in the velocity. Only if this is not the 
case, would the experimental evidence imply that fabric 
anisotropy exerts independent influence, that it is in fact 
a state variable - in addition to, and independent from, 
others such as the stress, density, velocity, and the gran- 
ular temperature. Generally speaking, of course, fabric 
anisotropy could well behave like many other local or mi- 
croscopic variables, is different in different circumstances 
but, being a dependent quantity, would be given if the 
complete set of state variables is. 

Fig. [1] demonstrates that the velocity difference can 
indeed be explained by the density difference. Although 
there is some residual discrepancy between theory and 
observation, many possible reasons remain, making the 
case for fabric anisotropy less than cogent: The calcula- 
tion was carried out under simplifying assumptions, and 
the experiment contains deviation from ideal conditions. 
The employed theory, GSH (for granular solid hydrody- 
namics) [H, contains only two elastic coefficients, leaving 
the energy independent of the third strain invariant, see 
section IIIII On the experimental side, possibly neither 
the density nor the static stress is uniform, as assumed. 
The latter is mainly a result of the elliptical geometry, 
and the inhomogeneity in density should especially be 
dependent on preparation. 

Having adjusted the parameters such that the density 



difference from the two preparations suffices to explain 
the velocity difference in the compressional waves, Khi- 
das and Jia find that the discrepancy in the shear wave 
velocities remains, and depends on the preparation, see 
their Fig. 8. Accepting this as an evidence for the in- 
dependent influence of fabric anisotropy, however, pre- 
sumes complete confidence in the density dependence of 
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FIG. 1: Symbols represent measured wave velocities versus 
the applied stress 033, as reported by Khidas and Jia, for two 
different preparations, decompaction and rain, corresponding 
respectively to the packing ratio of (j> = 0.606 and <f> — 0.643. 
(The geometry is triaxial, the preferred direction along the 
vertical is denoted as 3.) More specifically, rhombuses are 
vertical compressional waves c v ; circles horizontal ones ch- 
The remaining three are shear waves, with varying direction 
and polarization, see [jj. The lines are theoretical results, ob- 
tained by employing GSH with density-dependent elastic co- 
efficients, but excluding any possible independent influences 
from preparations. In the employed approximation of two 
elastic coefficients, assuming uniaxiality, GSH yields the same 
velocity for all three shear waves - a justified approximation, 
as their difference is clearly small. 
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emt, the theory employed in [TJ to account for the wave 
velocities. GSH sports a more realistic density depen- 
dence and the same discrepancy is small. In addi- 
tion, the neglected elastic coefficients associated with the 
third strain invariant should also depend on the density, 
further reducing the discrepancy. 

We derive the elastic modes from GSH in section [Til 
and apply the results to the experimental circumstances 
of [TJ in section Hill 



II. ELASTIC MODES IN THE GSH THEORY 

GSH theory is a hydrodynamic theory meant to ac- 
count for a wide spectrum of granular behavior, see 043 
for derivation and explanation. Its variables consist of 
the density p, momentum pvi, the elastic strain uy, and 
two entropies, the true and the granular one, s,s g . To 
consider elastic waves, we only need the equations 



d t (pvi) + Vj (cry + pviVj) 
uilSij 



d t u i3 - (1 - a) Vij -\ 



n 



0, 



= o, 



with the stress given as cry = (1 — a) itij — Cg^ll^ij ~Vg v ij > 
Try = —dw/duij, and w denoting the elastic energy. The 
coefficients a and t -1 ,Tj~ account respectively for the 
typical granular phenomena of stress softening and strain 
relaxation that give rise to plasticity. These coefficients 
and the viscosities £ s are all functions of the granular 
temperature and vanish for T g — > 0. For small wave 
amplitudes, the granular temperature T g is also small, 
making all these effects negligible, and granular media 
essentially elastic. This is the limit we shall consider here, 
while postponing the account of wave dissipation from 
viscosity and plasticity to a forthcoming paper. Taking 
a,r _x ,Tf 1 ,J7 fl ,Cg = 0, and uy = \ (d l U ] ■+ djUi), we 
have 



pUi+M ijkl V jk Ui = Q, 



(1) 



a wave equation for the displacements Ui, with the stress- 
dependent stiffness tensor Myy = d 2 w / du^duj-i- 
We take the granular elastic energy to be 



(2) 



w = BVA(2A 2 /5 + u 2 J£) , 



with A = — uu, u 2 = itytt?., and the superscript 
denoting the traceless part of any tensor, eg. tt^ = 
Uij + SijA/3. The two elastic coefficients are B and £, 
where £ is a constant, and B density dependent, 



B = B 



(3) 



where <j> = p/po is the packing ratio, with 4 , C p 1 4'ip re ~ 
spectively the ratio of closest and loosest packing, and 
4>\ p given by = (110 cp + 9cAf p )/20. 



Inverting the stress-strain relation Try = 7ry(ttfcj), tak- 



ing P = ttu/ 3 and n 2 = tAtA, we express the elastic 



strain as 



2L = 1 + Vl-^/(2P 2 ) 
A = (LP/B)$, u s =Tr s /(2fj,), 
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1/3 



(4) 
(5) 

(6) 



For uniform stress and strain, in the principal axes 
of both the stress and strain [see the first of Eqs flBJ)], 
?ry = nu\5ij, the stiffness tensor has the simple form 

Mijkl — — (^~~g~^J ( a (ik)$ijfikl ~ fiikdjl — SuSjk) , 
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where there is no summation over bracketed subscripts. 
So, for instant 



Mi 122 = 



B (LP 



i V b 
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A plane wave ansatz ~ exp{uit + kx) to eq.((T]) results in 
the eigenvalue problem (Ify — AcTy) t/ z = 0, where 

i£« = (ft,) + /3(i))riini/2 - S a , 
c 2 =u; 2 /k 2 = -XB(LP/B)i/(pO, 

2 



^ LP v 2 3 



4LP 



and rii — h-Jk. Eigenvalues A yield the velocity of 
propagation c, and the corresponding eigenvectors U 
the polarization. The general solution is given by the 
eigenvalues Ai = — 1 and A 2 ,3 = \ [r\ =F y/C ~ 2), where 

V = TLi^ andC^Eli^- 

As long as ft ^ ft ft> an d the wave vector is not 
along one of the principal directions, the associated dis- 
placements are 



Ui 



^2^3 (ft - ft) 

"3^1 (ft - ft) 

v nin 2 (A - ft), 

'«i(CtVC(/3i 
n 2 (CTVC(/3i 



(8) 



ft) + ftft) N 
ft) + ftft) 

ft) + ftft), 



In the case of propagation along principal directions the 
tensor K is diagonal from the very beginning, and we 
have, eg., for n\ = 1, Ai = ft — 1 and A2.3 = — 1 with 
Ui = (1,0,0), U 2 = (0,1,0) and U 3 = (0,0,1). Thus, 
plane waves traveling along the principal directions of the 
static stress field are either pure longitudinal or transver- 
sal modes. 
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For an isotropic background stress we have tt? = 0, 
and ft = ft = ft, and the compressional wave velocity 
shows the typical F« Hertz-scaling, 



'3 4\ J_B_ JP_ 



(9) 



For a uniaxial state, with 3 as the preferred direction, we 
have ft = ft and Ai = -1, A 2 . 3 = \ (r\ =F \J~C, - 2), with 
r, = (ft - ft) n? + ft, C = (ft 2 - ft 2 ) n\ + ft 2 . Depend- 
ing on n\ alone, A (and therefore also the velocities c) 
possesses rotational symmetry. The compressional wave 
velocity along 3, or n\ — 0, is given by 



c LP) ^ PG y b 



(10) 



The lateral compressional velocity, n\ = 1, is 



v v LP) & PG y b 



(ii) 



which clearly does not depend on B, only on £, deter- 
mined as £ = 3/4 from the measured velocities of both 
preparations, rain and decompaction. Clearly, £ does not 
depend much on the density, as assumed. Then, taking 
(j>i p = 0.555 and ft p = 0.664, any data point will fix Bo 
as 5.1 GPa. The result is Fig[TJ In fact, data points from 
compressional waves were used to fix £ and Bq, so the 
curves of shear wave velocities are free of fit parameters. 



The degenerate velocity for shear waves is a conse- 
quence of the simple form given for the energy in Eq ([2]) . 
Generally speaking, there are three strain invariants: 
A = -Tru, u 2 s ee Trtt 2 - ± (Trw) 2 , and u IH = (Tra) 3 - 
3Tr uTyu 2 + 2Tr u 3 , of which only the first two enter the 
expression of Eq @. Including the third invariant, we 
did find a discrimination of shear wave velocities under 
static shear. [Note that although um is of third order, 
its energy contribution, say ~ u////-\/A, may well be of 
the same order as the terms already in Eq ([2]) .] We are at 
present looking for a suitable extension of Eq taking 
into account in addition its effect on the Coulomb yield. 



Only for vanishing shear, tt®,tt s ~ > 0, L — > 1, do both 
formulas reduce to Eq Otherwise, the \[P- Hertz- 
scaling is not valid. The two shear wave modes for each 
of the two directions all have the eigenvalue, A = —1, so 



(12) 



III. ANISOTROPY OF GRANULAR SOUND 




IV. SUMMARY 



These results are now applied to analyze the data of [l[ . 
From Eqs (|10lll[) , we first calculate, for the uniaxial case, 
the ratio of vertical and horizontal compressional veloci- 
ties, 



c-h 



LPQ 



LPQ 



(13) 



We have shown that GSH, employing only two indepen- 
dent elastic coefficients, Bq = 5.1 GPa, £ = 3/4, is capa- 
ble of reproducing the experimental date of [l[ fairly well. 
Various possible reasons for the remaining disagreement 
are given and weighted, all rather more tangible than 
fabric anisotropy. Moreover, the density dependence of 
B and the density independence of £ is further validated. 
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